Intense P.1 (Gamma) diversification followed by rapid Delta substitution in Southern Brazil: a SARS-CoV-2 genomic epidemiology study.

The analyses of genetic traits, dispersion patterns and phylogenomics are essential for understanding the evolutionary forces driving SARS-CoV-2 viruses in these three years of COVID-19 pandemics. Brazil is one of the most affected countries in the world and not sufficient genomic studies have been performed. The emergence of P.1 lineage led to one of the most serious public health crises on record. Our study presents the genomic sequencing and characterization of 412 samples from Rio Grande do Sul state, in the Brazilian Southern region, during Gamma and Delta epidemic waves, in 2021. Additionally, molecular evolution tests were performed to identify positively selected sites in Brazil between 2020 and 2022, as well as offer some evolutionary perspective about the maintenance of multiple spike mutations in Omicron lineages. Genomic epidemiology analysis has indicated an intense P.1 (Gamma) diversification followed by rapid Delta substitution in Southern Brazil.

Genomic Surveillance of SARS-CoV-2 in Healthcare Workers: A Critical Sentinel Group for Monitoring the SARS-CoV-2 Variant Shift.

SARS-CoV-2 genome surveillance is important for monitoring risk groups and health workers as well as data on new cases and mortality rate due to COVID-19. We characterized the circulation of SARS-CoV-2 variants from May 2021 to April 2022 in the state of Santa Catarina, southern Brazil, and evaluated the similarity between variants present in the population and healthcare workers (HCW). A total of 5291 sequenced genomes demonstrated the circulation of 55 strains and four variants of concern (Alpha, Delta, Gamma and Omicron-sublineages BA.1 and BA.2). The number of cases was relatively low in May 2021, but the number of deaths was higher with the Gamma variant. There was a significant increase in both numbers between December 2021 and February 2022, peaking in mid-January 2022, when the Omicron variant dominated. After May 2021, two distinct variant groups (Delta and Omicron) were observed, equally distributed among the five Santa Catarina mesoregions. Moreover, from November 2021 to February 2022, similar variant profiles between HCW and the general population were observed, and a quicker shift from Delta to Omicron in HCW than in the general population. This demonstrates the importance of HCW as a sentinel group for monitoring disease trends in the general population.

Genomic characterization and molecular evolution of SARS-CoV-2 in Rio Grande do Sul State, Brazil.

SARS-CoV-2 is the virus responsible for the COVID-19 and has afflicted the world since the end of 2019. Different lineages have been discovered and the Gamma lineage, which started the second wave of infections, was first described in Brazil, one of the most affected countries by pandemic. Therefore, this study analyzed SARS-CoV-2 sequenced genomes from Esteio city in Rio Grande do Sul, Southern Brazil. We also comparatively analyzed genomes of the two first years of the pandemic from Rio Grande do Sul state for understanding their genomic and evolutionary patterns. The phylogenomic analysis showed monophyletic groups for Alpha, Gamma, Delta and Omicron, as well as for other circulating lineages in the state. Molecular evolutionary analysis identified several sites under adaptive selection in membrane and nucleocapsid proteins which could be related to a prevalent stabilizing effect on membrane protein structure, as well as majoritarily destabilizing effects on C-terminal nucleocapsid domain.

Reduced gut microbiota diversity in patients with congenital generalized lipodystrophy.

BACKGROUND: Previous studies suggest intestinal dysbiosis is associated with metabolic diseases. However, the causal relationship between them is not fully elucidated. Gut microbiota evaluation of patients with congenital generalized lipodystrophy (CGL), a disease characterized by the absence of subcutaneous adipose tissue, insulin resistance, and diabetes since the first years of life, could provide insights into these relationships. METHODS: A cross-sectional study was conducted with patients with CGL (n = 17) and healthy individuals (n = 17). The gut microbiome study was performed by sequencing the 16S rRNA gene through High-Throughput Sequencing (BiomeHub Biotechnologies, Brazil). RESULTS: The median age was 20.0 years old, and 64.7% were female. There was no difference between groups in pubertal stage, BMI, ethnicity, origin (rural or urban), delivery, breastfeeding, caloric intake, macronutrient, or fiber consumption. Lipodystrophic patients presented a lower alpha diversity (Richness index: 54.0 versus 67.5; p = 0.008). No differences were observed in the diversity parameters when analyzing the presence of diabetes, its complications, or the CGL subtype. CONCLUSION: In this study, we demonstrate for the first time a reduced gut microbiota diversity in individuals with CGL. Dysbiosis was present despite dietary treatment and was also observed in young patients. Our findings allow us to speculate that the loss of intestinal microbiota diversity may be due to metabolic abnormalities present since the first years of life in CGL. Longitudinal studies are needed to confirm these findings, clarifying the possible causal link between dysbiosis and insulin resistance in humans.

Detection of Intestinal Dysbiosis in Post-COVID-19 Patients One to Eight Months after Acute Disease Resolution.

The intestinal microbiota plays an important role in the immune response against viral infections, modulating both innate and adaptive immune responses. The cytokine storm is associated with COVID-19 severity, and the patient's immune status is influenced by the intestinal microbiota in a gut-lung bidirectional interaction. In this study, we evaluate the intestinal microbiota of Brazilian patients in different post-COVID-19 periods, and correlate this with clinical data and the antibiotic therapy used during the acute phase. DNA extracted from stool samples was sequenced and total anti-SARS-CoV-2 antibodies and C-reactive protein were quantified. Compared with controls, there were significant differences in the microbiota diversity in post-COVID-19 patients, suggesting an intestinal dysbiosis even several months after acute disease resolution. Additionally, we detected some genera possibly associated with the post-COVID-19 dysbiosis, including Desulfovibrio, Haemophillus, Dialister, and Prevotella, in addition to decreased beneficial microbes, associated with antibiotic-induced dysbiosis, such as Bifidobacterium and Akkermansia. Therefore, our hypothesis is that dysbiosis and the indiscriminate use of antibiotics during the pandemic may be associated with post-COVID-19 clinical manifestations. In our study, 39% (n = 58) of patients reported symptoms, including fatigue, dyspnea, myalgia, alopecia, anxiety, memory loss, and depression. These data suggest that microbiota modulation may represent a target for recovery from acute COVID-19 and a therapeutic approach for post-COVID-19 sequelae.

Emergence of Two Distinct SARS-CoV-2 Gamma Variants and the Rapid Spread of P.1-like-II SARS-CoV-2 during the Second Wave of COVID-19 in Santa Catarina, Southern Brazil.

The western mesoregion of the state of Santa Catarina (SC), Southern Brazil, was heavily affected as a whole by the COVID-19 pandemic in early 2021. This study aimed to evaluate the dynamics of the SARS-CoV-2 virus spreading patterns in the SC state from March 2020 to April 2021 using genomic surveillance. During this period, there were 23 distinct variants, including Beta and Gamma, among which the Gamma and related lineages were predominant in the second pandemic wave within SC. A regionalization of P.1-like-II in the Western SC region was observed, concomitant to the increase in cases, mortality, and the case fatality rate (CFR) index. This is the first evidence of the regionalization of the SARS-CoV-2 transmission in SC and it highlights the importance of tracking the variants, dispersion, and impact of SARS-CoV-2 on the public health systems.

Healthcare-Associated Infections-Related Bacteriome and Antimicrobial Resistance Profiling: Assessing Contamination Hotspots in a Developing Country Public Hospital.

Hospital-built environment colonization by healthcare-associated infections-related bacteria (HAIrB) and the interaction with their occupants have been studied to support more effective tools for HAI control. To investigate HAIrB dynamics and antimicrobial resistance (AMR) profile we carried out a 6-month surveillance program in a developing country public hospital, targeting patients, hospital environment, and healthcare workers, using culture-dependent and culture-independent 16S rRNA gene sequencing methods. The bacterial abundance in both approaches shows that the HAIrB group has important representativeness, with the taxa Enterobacteriaceae, Pseudomonas, Staphylococcus, E. coli, and A. baumannii widely dispersed and abundant over the time at the five different hospital units included in the survey. We observed a high abundance of HAIrB in the patient rectum, hands, and nasal sites. In the healthcare workers, the HAIrB distribution was similar for the hands, protective clothing, and mobile phones. In the hospital environment, the healthcare workers resting areas, bathrooms, and bed equipment presented a wide distribution of HAIrB and AMR, being classified as contamination hotspots. AMR is highest in patients, followed by the environment and healthcare workers. The most frequently detected beta-lactamases genes were, bla SHV-like, bla OXA- 23 -like, bla OXA- 51 -like, bla KPC-like, bla CTX-M- 1, bla CTX-M- 8, and bla CTX-M- 9 groups. Our results demonstrate that there is a wide spread of antimicrobial resistance due to HAIrB in the hospital environment, circulating among patients and healthcare workers. The contamination hotspots identified proved to be constant over time. In the fight for patient safety, these findings can reorient practices and help to set up new guidelines for HAI control.

Equivolumetric Protocol Generates Library Sizes Proportional to Total Microbial Load in 16S Amplicon Sequencing.

High-throughput sequencing of 16S rRNA amplicon has been extensively employed to perform microbiome characterization worldwide. As a culture-independent methodology, it has allowed high-level profiling of sample bacterial composition directly from samples. However, most studies are limited to information regarding relative bacterial abundances (sample proportions), ignoring scenarios in which sample microbe biomass can vary widely. Here, we use an equivolumetric protocol for 16S rRNA amplicon library preparation capable of generating Illumina sequencing data responsive to input DNA, recovering proportionality between observed read counts and absolute bacterial abundances within each sample. Under specified conditions, we show that the estimation of colony-forming units (CFU), the most common unit of bacterial abundance in classical microbiology, is challenged mostly by resolution and taxon-to-taxon variation. We propose Bayesian cumulative probability models to address such issues. Our results indicate that predictive errors vary consistently below one order of magnitude for total microbial load and abundance of observed bacteria. We also demonstrate our approach has the potential to generalize to previously unseen bacteria, but predictive performance is hampered by specific taxa of uncommon profile. Finally, it remains clear that high-throughput sequencing data are not inherently restricted to sample proportions only, and such technologies bear the potential to meet the working scales of traditional microbiology.

The presence of SARS-CoV-2 RNA in human sewage in Santa Catarina, Brazil, November 2019.

Human sewage from Florianopolis (Santa Catarina, Brazil) was analyzed for severe acute respiratory syndrome coronavirus-2 (SARS-CoV2) from October 2019 until March 2020. Twenty five ml of sewage samples were clarified and viruses concentrated using a glycine buffer method coupled with polyethylene glycol precipitation, and viral RNA extracted using a commercial kit. SARS-CoV-2 RNA was detected by RT-qPCR using oligonucleotides targeting N1, S and two RdRp regions. The results of all positive samples were further confirmed by a different RT-qPCR system in an independent laboratory. S and RdRp amplicons were sequenced to confirm identity with SARS-CoV-2. Genome sequencing was performed using two strategies; a sequence-independent single-primer amplification (SISPA) approach, and by direct metagenomics using Illumina's NGS. SARS-CoV-2 RNA was detected on 27th November 2019 (5.49 ± 0.02 log10 SARS-CoV-2 genome copies (GC) L-1), detection being confirmed by an independent laboratory and genome sequencing analysis. The samples in the subsequent three events were positive by all RT-qPCR assays; these positive results were also confirmed by an independent laboratory. The average load was 5.83 ± 0.12 log10 SARS-CoV-2 GC L-1, ranging from 5.49 ± 0.02 log10 GC L-1 (27th November 2019) to 6.68 ± 0.02 log10 GC L-1 (4th March 2020). Our findings demonstrate that SARS-CoV-2 was likely circulating undetected in the community in Brazil since November 2019, earlier than the first reported case in the Americas (21st January 2020).

Swab pooling: A new method for large-scale RT-qPCR screening of SARS-CoV-2 avoiding sample dilution.

To minimize sample dilution effect on SARS-CoV-2 pool testing, we assessed analytical and diagnostic performance of a new methodology, namely swab pooling. In this method, swabs are pooled at the time of collection, as opposed to pooling of equal volumes from individually collected samples. Paired analysis of pooled and individual samples from 613 patients revealed 94 positive individuals. Having individual testing as reference, no false-positives or false-negatives were observed for swab pooling. In additional 18,922 patients screened with swab pooling (1,344 pools), mean Cq differences between individual and pool samples ranged from 0.1 (Cr.I. -0.98 to 1.17) to 2.09 (Cr.I. 1.24 to 2.94). Overall, 19,535 asymptomatic patients were screened using 4,400 RT-qPCR assays. This corresponds to an increase of 4.4 times in laboratory capacity and a reduction of 77% in required tests. Therefore, swab pooling represents a major alternative for reliable and large-scale screening of SARS-CoV-2 in low prevalence populations.

How new molecular tools can help bugbusters: a Burkholderia cepacia complex outbreak investigation.

An outbreak of bloodstream infection (BSI) caused by members of the Burkholderia cepacia complex (Bcc) took place from March 2012 until April 2014 involving thirteen patients. AIM: To describe an outbreak investigation of BSI Bcc and showing how genetic sequencing tools contributed to confirm the hypothesis of extrinsic contamination proposed by an observational study. METHODS: The Infection Control Department revised and reinforced good practices of infusion therapy and catheter care, visits to affected wards, a case control study, and environmental screening based on the case-control findings. RESULTS: Data from the case-control study found an association of cases with central venous catheter (OR 1.36; CI 1.15-1.67) and intravenous cisatracurium use (OR 10.75; CI 1.67-68.89). Visits to the operatory block revealed problems related to the cold chain used for the preservation of thermolabile cisatracurium. We could not retrieve Bcc from environmental samples using classic microbiology. New samples from the same surfaces were obtained for genetic sequencing. Bcc was identified in the cooler box, refrigerator and reusable ice packages. CONCLUSION: Environmental screening using genetic sequencing proved to be a useful tool for confirming our hypothesis of extrinsic contamination raised by the case-control study.

One year cross-sectional study in adult and neonatal intensive care units reveals the bacterial and antimicrobial resistance genes profiles in patients and hospital surfaces.

Several studies have shown the ubiquitous presence of bacteria in hospital surfaces, staff, and patients. Frequently, these bacteria are related to HAI (healthcare-associated infections) and carry antimicrobial resistance (AMR). These HAI-related bacteria contribute to a major public health issue by increasing patient morbidity and mortality during or after hospital stay. Bacterial high-throughput amplicon gene sequencing along with identification of AMR genes, as well as whole genome sequencing (WGS), are biotechnological tools that allow multiple-sample screening for a diversity of bacteria. In this paper, we used these methods to perform a one-year cross sectional profiling of bacteria and AMR genes in adult and neonatal intensive care units (ICU and NICU) in a Brazilian public, tertiary hospital. Our results showed high abundances of HAI-related bacteria such as S. epidermidis, S. aureus, K. pneumoniae, A. baumannii complex, E. coli, E. faecalis, and P. aeruginosa in patients and hospital surfaces. Most abundant AMR genes detected throughout ICU and NICU were mecA, blaCTX-M-1 group, blaSHV-like, and blaKPC-like. We found that NICU environment and patients were more widely contaminated with pathogenic bacteria than ICU. Patient samples, despite the higher bacterial load, have lower bacterial diversity than environmental samples in both units. Finally, we also identified contamination hotspots in the hospital environment showing constant frequencies of bacterial and AMR contamination throughout the year. Whole genome sequencing (WGS), 16S rRNA oligotypes, and AMR identification allowed a high-resolution characterization of the hospital microbiome profile.

Unusual RNA plant virus integration in the soybean genome leads to the production of small RNAs.

Horizontal gene transfer (HGT) is known to be a major force in genome evolution. The acquisition of genes from viruses by eukaryotic genomes is a well-studied example of HGT, including rare cases of non-retroviral RNA virus integration. The present study describes the integration of cucumber mosaic virus RNA-1 into soybean genome. After an initial metatranscriptomic analysis of small RNAs derived from soybean, the de novo assembly resulted a 3029-nt contig homologous to RNA-1. The integration of this sequence in the soybean genome was confirmed by DNA deep sequencing. The locus where the integration occurred harbors the full RNA-1 sequence followed by the partial sequence of an endogenous mRNA and another sequence of RNA-1 as an inverted repeat and allowing the formation of a hairpin structure. This region recombined into a retrotransposon located inside an exon of a soybean gene. The nucleotide similarity of the integrated sequence compared to other Cucumber mosaic virus sequences indicates that the integration event occurred recently. We described a rare event of non-retroviral RNA virus integration in soybean that leads to the production of a double-stranded RNA in a similar fashion to virus resistance RNAi plants.

ASR5 is involved in the regulation of miRNA expression in rice.

KEY MESSAGE: The work describes an ASR knockdown transcriptomic analysis by deep sequencing of rice root seedlings and the transactivation of ASR cis-acting elements in the upstream region of a MIR gene. MicroRNAs are key regulators of gene expression that guide post-transcriptional control of plant development and responses to environmental stresses. ASR (ABA, Stress and Ripening) proteins are plant-specific transcription factors with key roles in different biological processes. In rice, ASR proteins have been suggested to participate in the regulation of stress response genes. This work describes the transcriptomic analysis by deep sequencing two libraries, comparing miRNA abundance from the roots of transgenic ASR5 knockdown rice seedlings with that of the roots of wild-type non-transformed rice seedlings. Members of 59 miRNA families were detected, and 276 mature miRNAs were identified. Our analysis detected 112 miRNAs that were differentially expressed between the two libraries. A predicted inverse correlation between miR167abc and its target gene (LOC_Os07g29820) was confirmed using RT-qPCR. Protoplast transactivation assays showed that ASR5 is able to recognize binding sites upstream of the MIR167a gene and drive its expression in vivo. Together, our data establish a comparative study of miRNAome profiles and is the first study to suggest the involvement of ASR proteins in miRNA gene regulation.

The source of the river as a nursery for microbial diversity.

Bacteria are highly diverse and ubiquitous organisms that play a key role as drivers for ecosystem processes. The application of NGS (next-generation sequencing technologies) for 16S analysis has been broadly used for understanding bacterioplankton composition and structure. Most of studies conducted on aquatic ecosystems with 16S NGS have been in seawater and lakes. A few studies using NGS have been conducted in river environments and have suggested the presence of a bacterial seed-bank. We performed 16S highly variable V4 region high-throughput analysis in the Sinos River, which is located in one of most important Brazilian industrial centers. This region has several contrasts in its environmental characteristics, presenting a longitudinal gradient of eutrophication and making it a remarkable study site for observing the dynamics of bacterioplankton. We demonstrated consistent evidence for the existence of a bacterial seed-bank and its longitudinal persistence. Seasonal shifts reinforce the importance of the source of the river in maintaining the bacterial seed-bank that spreads throughout the river. Therefore, the preservation of the source of the river is important not only for hydrologic reasons but also to maintain the microbial composition and the ecological integrity of the river.

The Wall-associated Kinase gene family in rice genomes.

The environment is a dynamic system in which life forms adapt. Wall-Associated Kinases (WAK) are a subfamily of receptor-like kinases associated with the cell wall. These genes have been suggested as sensors of the extracellular environment and triggers of intracellular signals. They belong to the ePK superfamily with or without a conserved arginine before the catalytic subdomain VIB, which characterizes RD and non-RD WAKs. WAK is a large subfamily in rice. We performed an extensive comparison of WAK genes from A. thaliana (AtWAK), O. sativa japonica and indica subspecies (OsWAK). Phylogenetic studies and WAK domain characterization allowed for the identification of two distinct groups of WAK genes in Arabidopsis and rice. One group corresponds to a cluster containing only OsWAKs that most likely expanded after the monocot-dicot separation, which evolved into a non-RD kinase class. The other group comprises classical RD-kinases with both AtWAK and OsWAK representatives. Clusterization analysis using extracellular and kinase domains demonstrated putative functional redundancy for some genes, but also highlighted genes that could recognize similar extracellular stimuli and activate different cascades. The gene expression pattern of WAKs in response to cold suggests differences in the regulation of the OsWAK genes in the indica and japonica subspecies. Our results also confirm the hypothesis of functional diversification between A. thaliana and O. sativa WAK genes. Furthermore, we propose that plant WAKs constitute two evolutionarily related but independent subfamilies: WAK-RD and WAK-nonRD. Recognition of this structural division will further provide insights to understanding WAK functions and regulations.

New insights into aluminum tolerance in rice: the ASR5 protein binds the STAR1 promoter and other aluminum-responsive genes.

Aluminum (Al) toxicity in plants is one of the primary constraints in crop production. Al³âº, the most toxic form of Al, is released into soil under acidic conditions and causes extensive damage to plants, especially in the roots. In rice, Al tolerance requires the ASR5 gene, but the molecular function of ASR5 has remained unknown. Here, we perform genome-wide analyses to identify ASR5-dependent Al-responsive genes in rice. Based on ASR5_RNAi silencing in plants, a global transcriptome analysis identified a total of 961 genes that were responsive to Al treatment in wild-type rice roots. Of these genes, 909 did not respond to Al in the ASR5_RNAi plants, indicating a central role for ASR5 in Al-responsive gene expression. Under normal conditions, without Al treatment, the ASR5_RNAi plants expressed 1.756 genes differentially compared to the wild-type plants, and 446 of these genes responded to Al treatment in the wild-type plants. Chromatin immunoprecipitation followed by deep sequencing identified 104 putative target genes that were directly regulated by ASR5 binding to their promoters, including the STAR1 gene, which encodes an ABC transporter required for Al tolerance. Motif analysis of the binding peak sequences revealed the binding motif for ASR5, which was confirmed via in vitro DNA-binding assays using the STAR1 promoter. These results demonstrate that ASR5 acts as a key transcription factor that is essential for Al-responsive gene expression and Al tolerance in rice.

isomiRID: a framework to identify microRNA isoforms.

SUMMARY: MicroRNAs (miRNAs) have been extensively studied owing to their important regulatory roles in genic expression. An increasingly number of reports are performing extensive data mining in small RNA sequencing libraries to detect miRNAs isoforms and also 5' and 3' post-transcriptional nucleotide additions, as well as edited miRNAs sequences. A ready to use pipeline, isomiRID, was developed to standardize and automatize the search for miRNAs isoforms in high-throughput small RNA sequencing libraries. AVAILABILITY: isomiRID is a command line Python script available at http://www.ufrgs.br/RNAi/isomiRID/.

Heavy metal-associated isoprenylated plant protein (HIPP): characterization of a family of proteins exclusive to plants.

Metallochaperones are key proteins for the safe transport of metallic ions inside the cell. HIPPs (heavy metal-associated isoprenylated plant proteins) are metallochaperones that contain a metal binding domain (HMA) and a C-terminal isoprenylation motif. In this study, we provide evidence that proteins of this family are found only in vascular plants and may be separated into five distinct clusters. HIPPs may be involved in (a) heavy metal homeostasis and detoxification mechanisms, especially those involved in cadmium tolerance, (b) transcriptional responses to cold and drought, and (c) plant-pathogen interactions. In particular, our results show that the rice (Oryza sativa) HIPP OsHIPP41 gene is highly expressed in response to cold and drought stresses, and its product is localized in the cytosol and the nucleus. The results suggest that HIPPs play an important role in the development of vascular plants and in plant responses to environmental changes.